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Efforts are underway to improve the design and sensitivity of gravitational waves detectors, with 
the hope that the next generation of these detectors will observe a gravitational wave signal. Such a 
signal will not only provide information on dynamics in the strong gravity regime that characterizes 
potential sources of gravitational waves, but will also serve as a decisive test for alternative theories 
of gravitation that are consistent with all other current experimental observations. We study the 
linearized theory of the tensor- vector-scalar theory of gravity (TeVeS) with generalized vector action, 
an alternative theory of gravitation designed to explain the apparent deficit of visible matter in 
galaxies and clusters of galaxies without postulating yet undetected dark matter. We find the 
polarization states and propagation speeds for gravitational waves in vacuum, and show that in 
addition to the usual transverse-traceless propagation modes, there are two more transverse modes 
and two trace modes. Additionally, the propagation speeds are different from c. 



I. INTRODUCTION 



Gravitational wave (GW) science is raising in popularity, with several GW detectors in operation in the United 
States (LIGO), Europe (VIRGO and GEO) and Japan (TAMA), and efforts under way to improve their design and 
O \ increase their sensitivity LISA, the laser interferometer space antenna, is expected to fly in the next decade, with 
LISA pathfinder marking the way Q • The detection of gravitational waves will convey new information on dynamics 
5^ ' of systems in the strong- field gravity limit, as well as in regions of the universe which are opaque to electromagnetic 
radiation. Gravitational waves will also serve as a test for alternative theories of gravity. Despite the outstanding 
success of General Relativity (GR) in the solar system, it fails to explain dynamics on galaxy and galaxy cluster 
scales without postulating large amounts of yet-undetected Dark Matter (DM). An alternative theory of gravitation 
J> ■ that does not require DM is TeVeS, the tensor- vector-scalar theory of gravity, which was designed as a relativistic 
lO ! implementation of Milgrom's Modified Newtonian Dynamics paradigm (MOND). MOND explains the asymptotical 
■ flatness of galaxy rotation curves without postulating the existence of yet undetected dark matter, as well as the 
sharpness of the Tully-Fisher relation which correlates luminosity of a disk galaxy with its asymptotic rotational 
velocity in a natural context. Milgrom proposed that Newtonian gravity progressively fails as accelerations drop 
below a characteristic scale Oo — 10~ 10 m/s 2 which is typical of galaxy outskirts, and assumes that for accelerations 
, of order Oq or well below it, the Newtonian relation a = — V<I>Ar is replaced by 



/i(|a|/oo)a=-V$jv, (1) 



where the function fl(x) smoothly interpolates between p,(x) = x at x <C 1 and the Newtonian expectation p.(x) = 1 
at x 3> 1. This relation with a suitable standard choice of p,(x) in the intermediate range has proved successful not 
^ ■ only in justifying the flatness of galaxy rotation curves in regions where acceleration scales are much below Oo, but 
also in explaining detailed shapes of rotation curves in the inner parts in terms of the directly seen mass, and in giving 
a precise account of the observed Tully-Fisher law , L cx V c 4 , predicting the relation L = (Ga T) _1 V c 4 . This sharp 
relation, while obtained naturally in the framework of MOND, requires quite a fine tuning of dark halo parameters 
to be explained by the dark matter paradigm [H . 

However, MOND alone is only a phenomenological prescription that does not fulfill the usual conservation laws, 
nor does it make clear if the departure from Newtonian physics is in the gravity or in the inertia side of the equation 
F = ma. Moreover, it is non relativistic, and as such it does not teach us how to handle gravitational lensing or 
cosmology in the weak acceleration regimes. To address these issues, Bekenstein designed the Tensor- Vector-Scalar 
theory of gravity (TeVeS) Q, a covariant field theory of gravity which has MOND as its low velocity, weak acceleration 
limit, while its nonrclativistic strong acceleration limit is Newtonian and its relativistic limit is GR. TeVeS sports two 
metrics, the "physical" metric on which all matter fields propagate, and the Einstein metric which interacts with the 
additional fields in the theory: a timelikc dynamical vector field, A a , and a scalar field, (j>. The theory also involves 
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a free function F, a length scale I, and two positive dimcnsionless constants k and K. The scalar field in TeVeS 
provides the additional gravitational potential for matter, whereas the vector field provides the desired light bending 
properties, in a fashion similar to the constant unit vector in Sanders' stratified theory Q. The TeVeS equations arc 
derived from an action principle, thus ensuring conservation of energy and momentum. 

One aspect of gravitational waves in TeVeS has been investigated by Sotani, who calculated oscillation spectra 
of neutron stars in the theory generated by perturbations of the fluid [6(, as well as by perturbations of the tensor 
and scalar fields Q, and compared them to the oscillation frequencies predicted in GR. Recently, Sagi @ calculated 
the PPN parameters for TeVeS and showed that there is a link between the cosmological value of the scalar field, 
4> c , and the coupling constant of the vector field, K, which prevents the scalar field from evolving with cosmological 
time as predicted by its equation of motion. This adds to existing evidence of dynamical problems in TeVeS with 
Maxwell-like vector action, that has been provided by Seifert Q and by Contaldi et al. [l(|. Consequently, in this 
work we investigate gravitational waves in vacuum for the generalized version of TeVeS introduced by Skordis , and 
determine the complete set of mode speeds and polarizations for generic values of the free parameters in the theory. 
This version of TeVeS has a vector action of the most general form which is quadratic in derivatives of the vector 
fields, and scalar and metric actions as in the original formulation of TeVeS, thus preserving the correct MOND and 
Newtonian limits. The vector field has four coupling constants associated with it, instead of just one in the original 
version of TeVeS. The results presented in this paper will be useful in determining the predictions of TeVeS with 
regard to gravitational wave emission from astrophysical sources. 

II. GENERALIZED TEVES EQUATIONS 

The metric, matter and scalar actions are given by Q 

S g = ^J g^R^yTg^x, (2) 
S m = j C (fa, f a ,f«,-- ') V^d*x, (3) 
Ss = ~2kk^G J ^(^h^M ^gd 4 x, (4) 

Above h Q/3 = g a @ — A a A@ with A a = g a ^Ap. In the scalar's action k is a dimensionless positive parameter while 
£ is a constant with the dimensions of length, and T a dimensionless free function. G is the gravitational coupling 
constant, and is not equal to the measured Newton's constant, which will be denoted Gn- The metric which couples 
to the gravitational fields is g^,, dubbed here the Einstein metric. The matter fields, which include all standard model 
particles, couple to the physical metric, which is composed of the three gravitational fields 

g a /3 = e~ 2ct >g al 3 - 2A a A p sinh(20). (5) 

The vector action is taken to be of the most general form quadratic in derivatives of the vector fields, as follows: 



S v = -j^J V~gd A x (-F a pF a P + -^S a0 S af3 + K 2 {VAf + K 4 A a A a - X (A a A a + 1)J (6) 

where F a p = A a -p — Ap ;a , S a p = A Q;( g + A^ a and A a = A^ A%. The Ki are dimensionless coupling constants. K is 
the coupling constant of the original version of TeVeS. 

Variation of the action with respect to g al3 yields the TeVeS Einstein equations for g a p 

G af} = 8nG (f al3 + (1 - e- 4 *) A»T^ a A p) + r Q/J ) + 6 a p, (7) 

where Vr a Ap\ = v a Ap + A a Vj3, etc. The sources here are the usual matter energy-momentum tensor T a p (related 
to the variational derivative of S m with respect to g a ^), as well as the energy-momentum tensors for the scalar and 
vector fields, 
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= kPlF S <j> -y(j) S . 



(9) 



(10) 



Each choice of the function J r (y) defines a separate TeVeS theory. Its derivative /i(y) functions somewhat like the \x 
function in MOND. For y > 0, fi(y) ~ 1 corresponds to the high acceleration, i.e., Newtonian, limit, while the limit 
< < 1 corresponds to the deep MOND regime. In the MOND regime, fi(y) ~ y/y/D, with D a dimensionlcss 
constant. We shall only consider functions such that T > and fj, > for either positive or negative arguments. 

The equations of motion for the vector and scalar fields are obtained by varying the action with respect to <f> and 
A a , respectively. We have 



[Ky)h a0 4>, a] . p = kG [g a ? + (1 + e" 4 *) A a A?] f a p , 

for the scalar and 

K\J a F al3 + K+\7 a S al3 + K 2 V P (V • A) - Kji'V^K + If 4 V„ (A^A'^j + AA& + ^-^A a (f>^g 
= 8nG (1 - e" 4 ^) g pa f ai A~< 
for the vector. Additionally, there is the normalization condition on the vector field 

A a A a = g a p A a AP = -l. 



(11) 



(12) 



(13) 

The A in Eq. ()12j) . the lagrange multiplier charged with the enforcement of the normalization condition, can be 
calculated from the vector equation. 



III. METRIC, VECTOR AND SCALAR PERTURBATIONS ON A CURVED BACKGROUND IN 

TEVES 

We start by considering metric, vector and scalar perturbations on a curved background. One can derive results 
on the structure of the solutions of the TeVeS equations, following the accounting system described in [l2|], Ch. 2. In 
order to discern the perturbation from the background we assume that in some coordinate system we can write the 
metric, vector and scalar as 

gap = 9aP + h a /3 (14) 

A a = A a + u a (15) 
= 0s + 6<t>, (16) 

where g a p, A a and <pB have a typical scale of variation Lb, on top of which small amplitude perturbations are 
superimposed, characterized by a scale l g satisfying £ g <C Lb (alternatively, the distinction can be made in frequency 
space, with the background characterized by a frequency much lower that the perturbation). Additionally, we assume 
that the background metric, vector and scalar are O(l), whereas the perturbations are of order e <SC 1. Writing the 
vacuum TeVeS Einstein equations in the form 

R a p = (SttGt^ + 6^) (S^p - \g a p 9 n, (17) 
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and expanding the Ricci tensor to 0(e 2 ), we can split the TeVeS equations into two parts: a low frequency, long 
wavelength part which describes how the background is affected by the perturbations, and a high frequency, short 
wavelength part which describes the propagation of the perturbations on the background, as follows 

\ Low 

(SnGr^ + 0„„) (8»8% - - 9a0 gn) (18) 

and 

/ i \ High 

Rap = -\R [ al] Hl9h + f (SttGt^ + 6^) {8*8} - -gapgn) (19) 

with Ofiu given by Eq. ([9]) and t^ v by Eq. l[8|). In the above, R a p is constructed from g a p, and contains only low 
frequency modes. is linear in h a p, and is thus high frequency, whereas R^p, which is quadratic in h a p, can 

contain both low-frequency modes generated from terms with nearly equal but opposite high wave-vectors, as well 
as high-frequency modes. In the absence of matter sources, the equations of motion force the amplitude of the 
perturbations e to be equal to £ g /Ls (this is true only in the absence of sources; if sources are present then they 
determine the background curvature whereas e <C £ 9 /Lb)- This comes about since the Ricci tensor contains terms 
which are quadratic in derivatives of the metric, so that on the left hand side of Eq. (|18[) we have terms quadratic 
in derivatives of the background metric, which are proportional to 1/L 2 B , and on the right hand side we have terms 
quadratic in derivatives of the perturbation, which are proportional to (e/£ g ) 2 . One can then set Lb = 1 and use 
a single expansion parameter \u a \, \S(j>\ ~ e ~ £ g -C 1. From Eq. (fT5| we see that terms quadratic in the 

perturbations induce changes in the background, whereas the propagation equations, Eq. (p~9|) . are first order in the 
perturbation. The leading term in the propagation equations, which contain second derivatives of the perturbations, 
is then 0(l/e), and first derivatives of the perturbations or of the background metric are 0(1). 

Regardless of the form of n(y), the scalar equation separates from the vector-tensor equations, since for all y 
< fi(y) < 1 and < J- < y. We can thus evaluate the order of magnitude of r a p by taking fi = 1 and T = y; this 
will give us a bound from above: r a p ~ </> jQ ,0^. Then 4> }a (f>,p ~ 0((e/£ g ) 2 ) = 0(1). The contribution of the scalar field 
to the vector equation, ^-fxA a (f> j0 ,gP' y (j)^, is of the same order of magnitude. Hence the scalar contributions to the 
equations for the vector and the metric are O(l), and are an order e smaller than the leading contributions, containing 
second derivatives of perturbations. Consequently, the vector-metric system separates from the scalar equation on a 
curved background, and one can treat the two separately. 



A. Vector-metric perturbations 



We are interested in the lowest order terms in the metric and vector equations, which are 0(l/e). To this order, we 
can approximate g a p by the Minkowski metric rj a p, and the background scalar and vector fields by A a w (1,0,0,0) 
and 4>b = 4>c ~ const. This approximation is valid only in the absence of matter, where the background curvature can 
be assumed to be almost absent. To avoid carrying factors of e 2 * c throughout the calculation, we chose coordinates 
in which the background Einstein metric is Minkowski; at the end of the calculation we will switch to coordinates 
in which the background physical metric is Minkowski. To first order in the perturbations, the difference between 
the two coordinate systems is only one of scale. For simplicity, we elected to work in the reference frame in which 
the vector field is at rest; when including matter content, the velocity of the matter frame with respect to the vector 
frame, v, would have to be accounted for, and the vector field would acquire a temporal component of the order of v 2 
and a spatial component of the order of v. We will work in units in which the speed of light is unity. 

We substitute Eqs. (|14|15ll6j) into the metric and vector TeVeS field equations in vacuum, Eq. (fl~7| and 

KV a F aP + K+V a S af) + K 2 V P (V • A) - K^V^A, + K 4 V [A p A a ^ + \A P + ^ M A Q 0, a /T0 )7 = 0. (20) 

To order 0(l/e), indices are raised and lowered with rjap, so that for example u l = Ui, etc. 

The temporal component of the vector equation gives A, which is first order in the perturbation, since the background 
lagrangc multiplier is zero. The spatial components are 

K 2 h jm + 2 A .//, ,.,„ + (K + K+- K 4 ) (h 00 ,oi - 2< o) + 2(# + K+K 3J = (21) 



5 



The metric equations are: 

00 :/;,,.,, - hu tjj = 2{K + K+- K A ) (h Qlfil - ^h 00 ,u + u] 0l ^j (22) 

01 :h Q: j d i + hjjjo - hoijj - hjjj = K 2 h jj . i0 + 2K + (hij tj0 + u\ J0 ) + 2(K 2 + K+)Uji (23) 

ij -hk(i,j)k - fro(j,i)o - h t ij - nhij - S i:j (/i o,oo + h k i M - UK) = K + (hij.oo + u j i0 + u l jQ ^ + K 2 Sij(hkk,oo + 2u* fc0 ) 

' (24) 

Round brackets denote symmetrization without a factor 1/2. h is the trace of the metric perturbation, and □ is the flat 
space d'Alembertian. The temporal component of the vector field perturbation is determined from the normalization 
condition: 

u° = h QQ (25) 

This system of equations is very similar to the system obtained in iEther linearized theory, therefore we will follow 
the analysis in [13( ■ 

As in GR, after choosing a frame within which the metric, vector and scalar fields have the form (fT4l [T5l [To]) , we 
are left with a residual gauge symmetry. Under an infinitesimal transformation of coordinates 

a? -W =x» + ? i {x), (26) 

the fields transform as 

K v =K v -{d^ v + d v ^), (27) 

=u" + Sb£ M , (28) 
54>' =6<f>. (29) 

If I^Ci^l ~ \h,j,u\; then the condition \h^ v \ 1 is preserved. We thus have the freedom to perform a linearized gauge 
transformation on the fields. The Lorentz gauge usually chosen in GR is of little use to us, since the additional terms 
in the equations that stem from the vector field stress-energy tensor are not simplified in this gauge. Instead, we 
choose as in [l3j to impose the four conditions hoi = and u\ = 0. It is easy to show that they can be obtained from 
a gauge transformation of the form (|2"6")l : starting from an arbitrary gauge that satisfies (Till' [in [To]) . one has to elect 
a vector £ M that satisfies 

fcoi = /HK-(&,o + eb,i) = (30) 

and 

< = Ui + &, = (31) 

Adding the vector equation and the spatial divergence of the metric equation, one gets a Poisson equation for £o; and 
can then be found by integrating the metric equation with respect to time. 

In the gauge in which hoi = and u\ = 0, the metric-vector system of equations takes the form: 

hij,ij — hujj = -{K + K + - K4,)hoo,u (32) 

hijjo — hjj.io = Kihjj^o + 2K+(hij t jo + (33) 

H(i,j)k - h ,ij - ah ij - % (Woo + hki,ki - Oh) = K + (h lh oo + u J i0 + u l j0 ^j + K 2 6ijhkk,oo (34) 

K 2 h jm + 2K+h ij , j0 + (K + K+ - K 4 ) {hoo,oi ~ 2t/ 00 ) + 2{K + K+)vi dj = (35) 

we now have thirteen equations and nine unknown functions (the divergence condition on the spatial part of the vector 
field leaves only two independent vector components), hence four of the equations above are redundant, and will serve 
as checks for our calculation. We will work with the six equations (f34| and with the three vector equations (|35p . 
since they are the equations controlling the dynamics of the system, whereas the rest of the equations are constraint 
equations. To solve the system, we assume plane wave solutions for the perturbations, in coordinates such that the 
wave vector is (ko, 0, 0, £3): 

hap = e a0 e ik » x " (36) 
u a = ce**"*". (37) 
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The equations become: 

u l : k k 3 K+e l3 + k 2 (K + K+- K A )e, t - k 2 3 (K + K+)e t = (38) 

u 3 : K 2 e u + (2K + + K 2 )e 33 + (K + K+ - if 4 )e o = (39) 

En : fc 3 2 eoo + ((1 + K 2 )kl - k 2 3 )e 22 + (2K+ + K 2 )k*e n + (1 + K 2 )k 2 e 33 = (40) 

E 22 : kjeoo + ((1 + K 2 )k 2 - fe|)en + (2K+ + K 2 )k 2 e 22 + (1 + K 2 )k 2 e 33 = (41) 

E 33 : (1 + K 2 )(en + e 22 ) + (2K+ + K 2 )e 33 = (42) 

E 12 : {{2K+ - l)k 2 + hi) ei2 = (43) 

E i3 : (2K+ - l)k 2 e l3 - 2K+k k 3 e l = (44) 

in all the above i = 1,2 and double indices imply summation. 

The above equations are homogeneous in and therefore there is no dispersion of the waves. If we define the 
wave speed to be s = ko/k 3 (it will be a real wave speed only if s 2 > 0), from the requirement that the determinant 
of this homogeneous linear system of equations is zero we obtain three wave speeds: 

1 



1 1 - 2K+ 

2 _ K + K+- 2KK+ 

82 ~ (1 - 2K+)(K + K+ - K 4 ) 
a = (K 2 + 2K + )(2-(K + K + ) + IU) 

' 3 (1-2K + )(K + K + -K 4 ){2 + 3K 2 + 2K+) [ ' 

Here we sec that for K + = 1/2 and for K + K + — K4 = the wave speeds are infinite. That K + K + — K4 = is 
disallowed is consistent with Skordis' result, that K — K 4 is the coefficient of the time derivative term in the vector 
cosmological perturbation equation fill ] , meaning that when K + K + — K4 = there is no growing mode in the vector 
field that can assist structure formation. For values of the coupling constants Ki for which s 2 is positive and finite, 
the dispersion relation obtained is linear, and |s| represents the propagation speed of gravitational disturbances. For 
these values, the theory has well defined propagating waves solutions. If s 2 is negative for a mode, then the frequency 
kg is imaginary, indicating the existence of exponentially growing or decaying solutions. In such a case the theory is 
unstable and hence presumably unphysical. 

The fields excitations corresponding to the wave speeds (|45|) arc: 



• Two transverse-traceless modes corresponding to s±: e 22 = —en 7^ 0, £12 7^ 0. 



Two transverse vector-tensor modes corresponding to s 2 : e % = \y — — 2K K+$i+K + -Ki) K+ ^ e ' 3 ' wn ^ n i = 1, 2. 



A trace mode involving the metric trace and the vector temporal component through the normalization condition, 

0(if 4 -(ifH 
K 2 +2K + 



j • j 1 \( T s , TS TS \ (1 + K 2 )(K4-(K + K + )) 

corresponding to s 3 : e = ^e 00 , en = e 2 2 = 2 {K + K+ - K A )t m , e 33 = ^ k 2 +2K+ e oo- 



The remaining metric equations agree with the above results. The modes are easily classified by their different 
propagation speeds, allowing us to naturally obtain the modes for the physical metric by simply substituting our 
result in the expression for the physical metric. To linear order, the physical metric is: 

5oo = e 2 ^(-l + h 00 -250) 
g Ql = -2u l sinh (2<f> c ) 

~ 9ij = e- 2 ^(<%(l - 28<t>) + hij). (46) 
Going to Minkowski coordinates through the following coordinate transformation 

x° = e^ x ° ia 3 =e -*. x i j (47) 

we get 

floo = -1 + h 00 - 26(f) 
g 0i = -2u l smh(2(j) c ) 

g ij = 8 ij {l-28(t>) + h ij . (48) 
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We must remember to transform the wave speeds to Minkowski coordinates as well: 



2 _ e- 4 ^ (K + K+- 2KK+) 
32 ~ (1 - 2K+)(K + K+ - K 4 ) 

2 _ e'^(K 2 + 2K+)(2 -JK + K+) + K A ) 
83 ~ (1 - 2K+)(K + K+- K 4 )(2 + 3K 2 + 2K+) 

We can then write the physical metric as f]^ u + , with the physical perturbation tensor given by 




cos [oj(t — z/si)] 






de 13 


de 23 













£13 


de 23 








£23 





£13 


£23 






[oj(t - z/s 3 )} 



cos [uj(t — zj ' s 2 )\ 



(50) 



(51) 



with d = 



(1-2K + )(2KK + -(K+K + )) 
K+(K+K+-K 4 ) 



sinh (2(j) c ), a=\(K + K+ - K 4 ) and b 



(1+K 2 )(K 4 -(K+K + )) 
K 2 +2K+ 



TcVcS exhibits the usual transversc-traccless propagation mode of GR, but at a speed different from c. Additionally, 
it has two more modes which excite both spatial and temporal directions, and one trace mode, which is not symmetric 
in all spatial directions. We will show in the next section that the scalar equation generates an additional trace mode. 

Incidentally, in the original version of TeVeS, which had all Ki but K zero, s\ = s 2 = e^ 2 ^" , showing no dependence 



on the coupling constant of the theor; 
This is consistent with the results in 



whereas s 3 = 0, meaning that the theory has no propagating scalar mode. 



B. Scalar perturbations 

Since the scalar field perturbations appear in the trace of the physical metric, we expect the scalar field to give rise 
to another, second trace mode. However, there is a subtlety in the derivation of the lowest order contribution of the 
scalar equation. The scalar equation in vacuum is 

[fx(y)h^ct>, a }. e = 0,. (52) 

To linear order in the perturbation, /j,(y) w /i(j/s) + A t '(ys)(y — Db)- With <f>B = 4>c ~ const., ub = and since for 
small y, fj,(y) ~ ^Jy, then /i'(ys) diverges. Therefore, to extract information on the propagation of scalar waves we 
have to relax the assumption 4>b ~ const, and allow it to depend on the space coordinates. This case has already 
been analyzed in using the WKB approximation, as is common practice for waves on a curved background; the 
scalar equation was shown to acquire the form 

= (h. af} + 2£,H a H f) ) S(j). a p (53) 

Ra s V'tB,, (54) 
^ v <t>B.p.4>B.M 

1 EE (55 ) 

amy 

here terms containing first derivatives of the scalar field perturbation were omitted, being O(l). Note that although 
originally this derivation was made assuming only scalar perturbations, and ignoring vector and metric perturbations, 
it is valid in the presence of metric and vector perturbations, owing to the separation of the metric-vector system from 
the scalar equation. Contributions of metric or vector perturbations to the scalar equation are also higher order, since 
there is no way of forming second-derivative terms of the metric or the vector perturbations in the scalar equation. 
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The wave speed was shown to be S4 < e 2 ^ B in the deep MOND regime, S4 = e 2< ^ >s /v / 2 in the Newtonian regime, 
and e^ B /y/2 < s 4 < y/1 + 2£e~ 2 ' i ' B /y/2 in the intermediate regime. Going back to the physical metric in Minkowski 
coordinates, Eq. (|48j) . we see that this generates an additional trace mode in the physical metric, at a speed different 
from that of the vector modes. In the notation of subsection (jlll A[) . it can be displayed as 



i- 


-25<j> 














-25<f) 








V 








-25(f) 














-26(/> 



cos [u)(t — z/s 4 )] 



(56) 



IV. CONCLUSIONS 



We investigated the propagation of gravitational perturbations in vacuum for the tensor-vector-scalar theory of 
gravity. We found that in the linear approximation on a curved background, the scalar equation separates from the 
vector-metric system of equations. We solved the vector-metric system of equations to lowest order in the background 
curvature, and obtained propagating wave solutions, with a linear dispersion relation, and three distinct wave speeds 
depending on the coupling constants of the theory and on the background value of the scalar field. The corresponding 
physical metric perturbations can be classified into a pair of transverse-traceless modes, another pair of excitations of 
the temporal-spatial components of the physical metric, and an asymmetric trace mode. Perturbations of the scalar 
equation were ill defined to lowest order in the background curvature; relaxing the restriction on the background 
scalar field, we were able to deduce from the analysis in Q that the scalar field gives rise to an additional trace mode, 
at a different propagation speed which depends on the background value of the scalar field and on the free function 
of the theory. 

TcVcS thus predicts six different modes of propagation for the gravitational field, at four distinct speeds, all different 
from the speed of light. Additionally, all speeds depend on a factor e _2 ^ B , which is expected to be close to unity when 
4>b = 4>c, but might induce a significant lag with respect to the speed of light for large values of the background scalar 
field. Such a lag might pose a problem for TeVeS; since the scalar and vector fields are coupled to matter via the 
physical metric, one would expect ultra high energy cosmic rays, whose velocity is close to c, to emit Cherenkov-like 
radiation of scalar and vector particles, if they move at a velocity higher than those scalar and vector particles. Such 
an emission would cause the cosmic rays to lose energy, and how much energy is lost would depend on the scalar and 
vector particles emission rate, on the distance traveled by the cosmic rays from their sources, and on the strength of 
the matter-fields coupling. Such Cherenkov-like radiation might pose very stringent restrictions on TeVeS' parameters, 
as was the case for Einstein- iEther theory (l4| . 

That problem could be avoided, at least for the physical perturbations originating in the vector-metric system, if the 
values of the coupling constants of the theory were such as to make the propagation speeds larger than c. The question 
of whether superluminal propagation in a theory with two metrics is allowed is still open jl5l - ll7l ] . If superluminal 
propagation could be allowed in TeVeS without disturbing causality, then one could also think of relinquishing the 
unconventional kinetic term in the action of the scalar field, which was introduced to prevent faster than light scalar 
waves. 
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